LOWER BOUNDS FOR THE TOTAL STOPPING TIME OF 3X + 1 

ITERATES 
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Abstract. The total stopping time aoo (n) of a positive integer n is the min- 
imal number of iterates of the 3a; + 1 function needed to reach the value 
1, and is +oo if no iterate of n reaches 1. It is shown that are infinitely 
many positive integers n having a finite total stopping time <Too(l), such that 
foo (n) > 6. 14316 log n. The proof involves a search of 3a; + 1 trees to depth 60, 
A heuristic argument suggests that for any constant 7 < -ysp ~ 41.677647, 
a search of all 3a; + 1 trees to sufficient depth could produce a proof that 
there are infinitely many n such that <Too(«) > 7logn. It would require a 
very large computation to search 3a; + 1 trees to a sufficient depth to produce 
a proof that the expected behavior of a "random" 3a; + 1 iterate, which is 
7 = log 4/3 ~ 6.95212, occurs infinitely often. 
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1. Introduction 

The 3a; + 1 problem concerns the behavior under iteration of the 3a; + 1 function 
T:Z^Z given by 



(1) T{n) 



§ if n = (mod 2) 
2^±i if n = 1 (mod 2) 



The notorious 3x+ 1 Conjecture asserts that for each n > 1, some iterate T^ k \n) = 
1. It has been verified for all n < 2.702 x 10 16 by Oliveira e Silva whose sub- 
sequent computations have extended this bound to 1.000 x 10 17 . The dynamical 
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behavior of the 3x + 1 function has been extensively studied. It is a determinis- 
tic process which nevertheless appears to exhibit a kind of pseudorandom behav- 
ior. Surveys on the 3a; + 1 problem can be found in Lagarias (|], Miiller jl0| and 
Wirsching [[L5| . 

To study the behavior of iterates of the 3a; + 1 function on the positive integers 
quantitatively, we make the following definitions. For each n we call the minimal 
k such that T^ k \n) — 1 the total stopping time of n and denote it tr m (n), letting 
o~oo(n) = co if it is otherwise undefined. A rescaled version of the total stopping 
time, which we call the stopping time ratio 7(77,), defined by 

2 7 (n) := -j . 

logn 

A trajectory converges if it has a finite total stopping time, and we consider it to 
consist of n and all its iterates up to the first iterate k with T^ k \n) = 1. Its parity 
sequence v(n) is the zero-one vector of length o~oo{n) giving the residue (mod 2) 
of n and its iterates up to T^-^n). The ones-ratio p(v) of a zero-one vector 
v is the ratio of the number of ones in the vector to its length. By extension 
the ones ratio p(n) of a convergent trajectory is defined to be the fraction of odd 
integers appearing in its parity sequence, i.e. the ones-ratio p(v(n)). For example, 
the trajectory of n = 3 is (3,5,8,4,2,1), with total stopping time <Too(3) = 5, 
parity sequence v(3) = (1,1,0,0,0) and ones-ratio p(3) = | = 0.4. In this case 
7(3) = Tsfa « 4.5512. 

This paper is concerned with rigorous results concerning the behavior of the 
stopping time ratio 7(n). It is well known that "f{n) > j^-j for all n, and that 
equality holds exactly for n = 2 k for k > 1. In this paper our object is to get 
lower bounds for the stopping time ratio, of the form 7(71) > 7, for some constant 
7 > , that hold for infinitely many n. We note that study of the size of the 
stopping time ratio is essentially the same as studying the size of the ones-ratio, as 
given by the following bounds. For any convergent trajectory one has 

(3) 7(n) > ' 



log 2 - p{n) log 3' 
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while if p(n) < .61 then, for any positive e, 



1 

7(n) " log2-p(n) log 3 + f : 

holds provided n is sufficiently large, with n > no(e). (It is believed that the 
condition p(n) < .61 should hold for all sufficiently large integers, see below, while 
to get a uniform bound in terms of e as above, one needs only assume p(n) < co < 
« .63092.) Our actual analysis will be based on study of the ones-ratio. 
As background, we review the current conjectures about the size of the stop- 
ping time ratio. These conjectures were suggested by theoretical results proved for 
stochastic models; these models make predictions about the total stopping time 
function which successfully match experimental evidence for iterates of the 3x + 1 
function. Simple models for average-case behavior of the stopping time function 
appear m Crandall §, La garias ||, Rawsthorne [Q, Lagarias and Weiss Q and 
Borovkov and Pfeifer ||. These models predict that the average value of 7(n) 
should be 

6.95212, 



log 4/3 



which corresponds to a ones-ratio of ^, and that nearly all values of 7(n) should 
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be within Q{{ * = of this value. This is well supported by experimental 



evidence. 

Lagarias and Weiss Q also developed stochastic models to describe the maximal 
values attained by the total stopping time. These models are are analyzed using the 
theory of large deviations. They predict that the limit superior of 7(71) as n — > 00 
should be a certain constant jbp ~ 41.677647, which corresponds to a ones-ratio of 
about 0.609091. (This constant is the solution to a certain functional equation given 
in ||.) The occurrence of values close to the extremal one are extremely rare events 
whose probability can be estimated in the stochastic models, see [|[ Section 6] ) . For 
example, the largest value to be expected for n < 10 17 is about 32, and Oliveira 
e Silva's record number is n = 1,008,932,249,296,231 with cr^n) = 1142 and 
7(n) « 33.0558. Backwards search methods have uncovered some larger integers 
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which have larger values of 7(71). V. Vyssotsky [jl3| found that 7(71) ~ 35.2789 for 

n = 37, 664, 971, 860, 959, 140, 595, 765, 286, 740, 059, 

which has croo(n) = 2565. He also found the current record value for 7(71) of about 
36.40, coming from a number around 10 110 . 

We now turn to the problem considered here, that of obtaining bounds a such 
that there are provably infinitely many positive n having 7(71) > 7. The current 
best value of a is that associated to the family of numbers n = 2 k — 1. It is well 
known that after k iterations one reaches T^{n) = 3 fe - 1. Since this number 
cannot get to 1 any faster than a power of 2 (which in fact occurs for k = 2), one 
obtains the rigorous lower bound 



(4) 7 (2 fc - 1) > ° ° « 3.729. 



log 2 + log 3 

This corresponds to a ones ratio about 0.387. It is expected that the numbers 3 fc — 1 



for most k have iterates behaving like a "random" integer, in which case the lower 
bound (Q) could be significantly improved for infinitely many k. However, as far as 
is currently known, it could be the case that for each e > the bound 

^oc(3 fc -l)<(7^+6)log(3 fc -l), 

holds for all sufficiently large fc; if so then (Q) could not be improved upon, asymp- 
totically as k — > 00.. 

The object of this paper is to improve this lower bound, as follows. 

Theorem 1.1. There exist infinitely many n > 1 having a convergent trajectory 
with a ones-ratio of at least || « 0.4827. Thus there is an infinite set of positive 
integers n such that 7(71) is finite and 

29 

(5) 7(n) > — — - — — — - w 6.14316. 

w n ' - 291og2 - 141og3 

This result is proved by a study of all 3a; + 1 trees to depth 60, by extensive 
computation, involving paths with up to 29 ones. In the rest of the paper we 
describe 3x + 1 trees, the form of proof certificates for lower bounds of the kind in 



Theorem 1.1 and the computations. 
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2. 3X + 1 Trees 

The 3a; + 1 iteration when run backwards from any fixed integer a produces a 
tree of preimages of a. 

A 3x + 1 tree Tk(a) is a rooted, labelled tree of depth k, representing the inverse 
iterates T~ 3 (a) for Q < j < k. The inverse map T (n) is multivalued: 

{{2n} if n = or 1 (mod 3), 

{2n,^-} if 7i = 2 (mod 3). 

The root node a is at depth 0, and a node labelled n at depth d of the tree is 
connected by an edge to a node labelled T{n) at depth d — 1 of the treefj. Thus 
the depth of a node is the number of edges in a path from it to the root node. As 
described in ||, the nodes labelled n = (mod 3) give rise only to a linear chain of 
nodes labelled n' = (mod 3) at higher levels. It is convenient to remove all such 
nodes and study a "pruned" tree T£{n) consisting of nodes n ^ (mod 3). That 
is, we study the inverse map on the set of positive integers n ^ (mod 3) given by 

{{2n} if 7i = 1,4,5 or 7 (mod 9), 

{2n, if 7i = 2 or 8 (mod 9). 

Figure [l] presents some examples of 7k (a) and T£{a). (Nodes n = 5 (mod 9) are 
circled to indicate that they have some preimage T _1 (ti) = (mod 3), and nodes 
7i = (mod 3) are indicated with a square.) 

We assign to each edge of a 3a; + 1 tree between depth d and d — 1 an edge label 
or 1, which is the value ti (mod 2) of the endpoint node at depth d. Traversing 
a path from a leaf node labelled ti in T^(a) at depth d to the root node a, the 
successive edge labels encode the initial part of the parity vector v(n) of n, and the 
zero-one weight vector w(ti; a) is the reversal of the vector v(n) to the root and 
gives the sequence of edge labels from the root node to the node n. The weight w 
of such a path is the sum of the edge labels on the path, and we sometimes say that 
the node n has weight w(n). Note that w ^ is the ones-ratio of this vector along 

1 We adopt a convention of "unrolling" any cycles under T, so that the same node label may 
appear at different levels of the tree if a cycle is present, cf. Figure hi 
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the path. The max-weight w max (a) of a tree T£{a) is the largest weight among all 
paths from a leaf node to the root node. Thus Wma *( a ) i s the maximal ones-ratio 
among any path in the tree T£(a) from a depth k node to the root. 

The structure of a pruned 3a; + 1 tree T^{a) of depth k is its isomorphism 
class as a rooted, edge-labelled tree. That is, we say that two pruned 3a; + 1 
trees T£{a) and T£(b) have the same structure if they are isomorphic as rooted 
trees by an isomorphism that preserves edge labels. The structure of a "pruned" 
3a; + 1 tree T£{a) is completely determined by knowledge of a (mod 3' +2 ), where 
I = w max (a). To demonstrate this, we grow the tree from the root, labelling each 
new node n (mod 3 m ), where m is determined recursively, as follows. The root 
node is assigned m = w max {a) + 2. For each node n (mod 3 m ) at level d, we assert 
m > 2, and if n = 1,4, 5 or 7 (mod 9) there is a single edge labelled to a node 
at level d + 1 labelled 2n (mod 3 m ). If n = 2 or 8 (mod 9), there is one node 
at depth d + 1 as above plus a second edge labelled 1 to a node at depth d + 1 
labelled (2n — l)/3 (mod 3 m_1 ). Since we started with m = w max {a) + 2 at the 
root node, and since no path has more than w max (a) edge labels of 1 , the value of 
the exponents m stays at least 2 all the way to, and including, the leaves at level k. 

We can improve on this bound, in one circumstance. We say that a tree T£(a) 
is critical at max-weight (or max-level) / if it has max-weight w ma x(a) — I and if 
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(i)T 5 (l) (ii)T 5 (4) (iii)T 5 *(4) 

Figure 1. 3x + 1 trees T k (a) and "pruned" 3a; + 1 tree T£{a) 
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all nodes at depth k — 1 have weights w max (a) — 1 or less in their path to the root. 
The structure of a critical tree is completely determined by a (mod 3' +1 ) where 
I = w max (a), rather than a (mod 3' +2 ). The proof is the same as before, noting 
that the exponent m > 2 on every node up to and including depth k—1, but m = 1 
may occur on some leaf nodes, on paths of weight exactly w max (a). 

For trees of depth k one has I < k, and the case I = k is necessarily critical, 
so by the discussion above all pruned tree structures 7jj,*(a) are determined by 
a (mod 3 fe+1 ). Thus there are at most 2 ■ 3 k distinct pruned tree structures T£{a) 
of depth k. The actual number R(k) of distinct tree structures is smaller but still 
grows exponentially in k. 

We now suppose given as data a congruence class a (mod 3' +1 ), and from it we 
may grow a "pruned" 3a; + 1 tree with root a to the (unique) depth k at which it 
becomes critical with max-weight /. We represent a (mod 3' +1 ) as a ternary (or 
3-adic) vector c = c(a) = (co, c\..., cj) of length I + 1, where 

i 

a = (qq_i...c ) 3 = ^Cj-3 J , 

3=0 

with Cj G {0, 1,2}. We label this tree 7^j(c). Using the procedure above, we can 
construct node labels for this tree consisting of such 3-adic vectors of variable length 
to at each node. Figure^ illustrates this, depicting the "pruned" tree T^(A), which 
is critical with max-weight 2, and the same tree 7j 2 ] (Oil) computed using root value 
4 (mod 27). 

We define the ones-ratio p(c) for a critical tree 7m (c) to be 



where k is the depth of the tree, and I is its criticality level. 

3. Lower Bound Certificates 

We describe a finite set of information ( "certificates" ) that are sufficient to yield 
proof that there are infinitely many n > 1 having a finite total stopping time whose 
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(ii) 73,(011) 



FIGURE 2. "Pruned" 3x + 1 tree 7[*j(c) with c = 110, for a = (011) 3 



ones-ratio satisfies the lower bound 



(6) 



p(n) > a . 



Definition 3.1. A prefix code C for ternary sequences is a finite set of codewords 
c = (c , Ci, . . . , q) G {0, 1, 2}' +1 of varying lengths I = 1(c) having the properties: 

(i) If c, c' G C with c^c', then c is not a prefix of c' or vice- versa. 

(ii) Any infinite sequence s = (s\, S2, S3, ... ) G {0, 1, 2} N ° has an initial segment 
(si, S2, . . • , sj+i) that is a codeword in C. 



Definition 3.2. A (lower bound) certificate for ones-ratio a consists of: 

(i) A finite list C of residue classes o (mod 3 i+1 ) for various lengths, I = 1(a), 
which are exhaustive, in the sense that their associated vectors c(a), together 
with the additional vector Co = 0, form a (ternary) prefix code. 

(ii) With each a G C with a (mod 3 /+1 ) is given a critical path vector w(a) G 
{0, l} k , where k satisfies 

l/k > a , 

such that w(a) has weight at most /, and if it has weight exactly I, then a one 
occurs in its last coordinate, and w(a) = w(n; a) for some leaf node of this 
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Table 1 . Certificate for p = | . 



3a; + 1 tree, i.e. w(a) is a labelled path in the (3a; + 1) tree of a (mod 3 !+1 ), 

from the root node to some node at depth k. 
The certificate vector certifies that the associated critical tree 7j^(c(a)) with max- 
level I has depttj^] k' with k' < k and has ones-ratio l/k' > a. 

As an example, we present in Table [l]a certificate for level a — 1/3. In this table 
the numbers a (mod 3' +1 ) are represented in their ternary expansion, so form a 
suffix code, rather than a prefix code. (The vectors c(a) reverse the order of the 
digits.) This certificate has max-weight 4 and max-depth 12, and consists of 12 
different trees. 

Definition 3.3. A strong lower bound certificate for ones-ratio a consists of a finite 
list C as above satisfying: 

(i') Same as (i) above. 

(ii') With each a (mod 3' +1 ) £ C are given two critical path vectors Wj(a) 6 
{0, l} kj , with weights lj < I, respectively, such that neither is a prefix of the 
other, and where kj and lj satisfy 

Ij/kj > a , 

2 A single path is a 3x + 1 tree sometimes can be grown to a depth greater than k, from 
knowledge of a (mod S'"*" 1 ), as long as it has weight at most I, and if exactly weight I then it has 
a one in its last coordinate. 
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and if has weight exactly I then it has a 1 in its last coordinate, and each 
v/j (a) can be found as a labeled path of length kj from the root in the 3x + 1 
tree of a (mod 3' +1 ). 

As an example, there exists a strong lower bound certificate for a — | which has 
max-weight 4 and max-depth 12, the same as the lower bound certificate above, 
but which consists of 42 trees. It appears in Table || below. 

Theorem 3.1. (i) The existence of a lower bound certificate for ratio a implies: 
Given any a > 1 with a ^ (mod 3), not in a cycle, there exist infinitely many n 
having 

(7) T^(n) = a, 

for some k > 1, and such that the ones-ratio p(n; a) := p(w(n; a)) of the path from 
n to a has 

(8) p(n; a) — l/k > a , 

(ii) The existence of a strong lower bound certificate for ones-ratio a implies: Given 
any a > 1 with a ^ (mod 3) the number of n < x with T^- k \n) — a for some 
k > 1 and with 

(9) p(n; a) — l/k > a , 

is at least c{a)x s for certain c(a) > and S — 5(a) > 0. 

Proof, (i) Since a is not in a cycle, all the nodes in the (infinite) "pruned" 3x + 1 
tree T*(n) with root node a have distinct values. 

We can find an infinite path s = (s\, S2, S3, . . . ) € {0, 1} N ° from the root node 
a in T*(a) which contains infinitely many nodes satisfying (^|). View a as a 3-adic 
integer, encoded as 

c(a) = (c , Cl ,c 2 ,c 3 ,...) e {0,1,2} N ° ; 

in which Co 7^ and Cj = for all sufficiently large j. Then c(a) has a prefix in the 
certificate C, of length I + 1, say. By property (ii) of the certificate, there exists a 
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Table 2. Strong certificate for p = | 



path in the associated tree 7^*](c(a)) from leaf node n to a with 

T( fel )(m) = a 



and 
(10) 



p(n;a) = -i- > a , 
ki 
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where k\ is the depth of 7j*j(c(a)). Now n\ £ Z>o is not in a cycle and rii ^ 
( mod 3), so we can repeat the same construction; finding a prefix of c(ni) of length 
l 2 + 1, say, with tree 7^ (c(ni)) of depth having a path from leaf node n 2 to ni 
with 

T^(n 2 ) =m 



and 



Thus 



p(n 2 ,ni) = ^->a. 
k 2 



T^ +k2 \n 2 ) =a 



and 

(11) p(^2, a) > min 



fa h 
k 2 ' fci 



> a 



Continuing inductively we find an infinite chain a *— n± <— *— TI3 *— . . . for 
which (0) and (§) hold. 

(ii) We use a similar argument. Suppose first that a is not in a cycle. Then 
at the first round we find two different paths in 7|*,(c(a)) with ni — > a, — > a 
satisfying (ph. At the j-th stage we produce 2- 5 distinct elements rij —> a satisfying 
(||). Furthermore, if the largest depth tree 7™ (c) for c G C is d then, all the elements 
rij are at depth at most jd from the root node a, hence each such element satisifes 
rij < 2i d a. Let 7r"(x) count the number of elements < x which satisfy (0) and (|J). 
It follows that for 2 jd a <x< 2^ +l 1 d a we have 

(12) <(x)>2^>i(^) 1/d = c( a )x 1 / d 

with c(a) = \a- x l d . 

Now suppose that a is in a cycle, i.e. is periodic. Then in the first tree 7j*j(c(a)) 
at least one of the two preimages n\ and n! 1 of a that satisfies (^|) cannot be in the 
periodic cycle containing a. Let it be n\, and we may then apply the argument 
above to the tree with root n\ for which 

K 1 ( x )>l(^r 1/dxl/d ■ 
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But any such element n — > n\ — > a by adjoining the path from n\ to a, and the 
inequality ( [Tlj) gives 

p(n,a) > mm[p(n,ni), p(ni,a)] > a . 

Thus 

(13) > > i^O-V^V^ > I -VV/ d , 

since ni < 2 d a. □ 



Theorem 3.1 reduces the problem of finding lower bounds for pin) for infinitely 
many n to that of finding suitable certificates, resp. strong certificates. Such 
certificates can be searched for by computer, as we describe next. We note that for 
any fixed a > \ the stochastic models studied in Lagarias and Weiss || predict that 
the number of n < x that have p(n) > a is bounded above by a; 7 ( Q ) +0 M for a certain 
exponent 7(a) < 1. (The value 7(a) is effectively computable as the fixed point of 
a certain functional equation.) It follows that strong lower bound certificates for 
such a, if they could be found, would produce a lower bound qualitatively of the 
correct form as what would be expected from the stochastic model predications, 
but (presumably) with a much smaller exponent than 7(a). 

4. Computational Results 

We may search for lower bound certificates resp. strong lower bound certificates 
by a "greedy" algorithm, as follows. 

Certificate Search Algorithm. 

(0) Input Data: a with < a < 1. 

(1) Initialization: C = {c^ : i G {1,2}, j £ {0,1,2}} with = £ {0,1, 2} 2 . 
Set level 1 = 1. Declare all vectors cy G C open. 

(2) Tree Test Step: Select an open vector c £ C, and let Z + 1 be the length of 
c. Compute T[*(c) and determine its depth d = d(l,c). If l/d > a declare 
the vector c closed. Otherwise split c into three vectors of length / + 2, (c, 0), 
(c, 1) and (c, 2), at level I + 1. Update C by deleting c and adding the three 
new vectors. Declare the three new vectors open. 
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(3) Termination Test: If all c G C are closed, halt. Otherwise, repeat the tree 
test step. 

In order to keep the set of open prefixes C as small as possible during this 
algorithm, it is convenient to choose at each test step to examine an open prefix 
vector that is one of those currently of greatest level. 

This algorithm is not guaranteed to halt, but if it does, it produces a lower 
bound certificate for ratio a. The correctness of the algorithm is based on the 
observation that C is a prefix code at all times: the "splitting" step replacing c by 
(c, 0), (c, 1) and (c, 2) preserves the prefix code property. It is also easy to see that 
the algorithm will find a certificate C if any certificate exists for ones-ratio a, and 
it will be one of minimal depth. 

The algorithm above is easily modified to give a Strong Certificate Search Algo- 
rithm. It is exactly the same, except that the rule for closing a node is modified: A 
vector c S C at level / is declared closed (only) if l/d > a and at least two paths of 
weight I exist from the root to a leaf of T*(c). If the algorithm halts, it produces 
a strong lower bound certificate for ones-ratio a. 

We searched by computer for certificates and strong certificates, for different 
values of a. The certificate and strong certificate for level a — 1/3 presented in 
Table [l] and Table || were found in this way. The size of certificates grows rapidly 
with increasing a. Table |§| presents data giving the maximal value of a for which a 
certificate exists in which all trees have a critical path containing at most I ones, for 
I < I < 29, as well as a (not necessarily maximal) value of a for which a certificate 
exists for Z = 29. 

Theorem 4.1. A lower bound certificate C exists for ones- ratio a = || « 0.4827. 

Proof. The certificate was found by computer search, of size indicated in Table [| 
The computer search required examining trees to depth 60, with some paths having 
29 ones, although the path with the largest ones ratio occured in a tree of depth 
58, having 28 ones. It is an interesting feature of this certificate that the "worst" 
tree is not one of maximal depth in the search. □ 
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Table 3. Size and depth of certificates 



Proof of Theorem 1.1. Since a = 1 is in a cycle, we cannot apply Theorem 3. 1 (i) 
directly. Instead we consider n — 41 and note that p(n) > and j(n) > 20. Now 



we apply Theorem 3.1 (i) to n = 41 using a = || and the certificate of Theorem 4.1 



Then the elements n — > 41 produced in Theorem 3.1 (i) have n — > 41 — > 1 hence 

14 

> min[p(n;41),p(41)] > — , 

as required. □ 
The corresponding results for strong certificates are given in Table ^, giving 
the maximal value of a for which a strong certificate exists having trees with all 
critical paths having at most I ones. This table is slightly less extensive because 
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Table 4. Size and depth of strong certificates 



larger searches were required. The breakpoint values of a are not the same as for 
certificates. In many cases they do not involve a larger depth for the largest tree, 
but sometimes do. 

Let 7r a (x) count the number of n with 1 < n < x which have finite total stopping 
time and ones-ratio p(n) > a. 



Theorem 4.2. A strong certificate C exists for ones-ratio a = g| w 0.4745. Fur- 
thermore there is a positive constant c such that 



ir a (x) > cx 1/59 . 
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Proof. The proof is similar to that for Theorem 4.1. The bound for the exponent 
2,1/59 comes 

from ©. □ 

We remark on the size of the certificates versus the computation needed to find 
them. Although the certificates are large, they can be verified without computing 
the entire tree T^(c); one simply backtracks along the given test path(s). However 
to find a minimal depth certificate, it appears necessary to calculate much of the 
tree 7^* (c) . Thus a much larger computation seems needed to find the certificate 
then to check it afterwards. Note that Lagarias and Weiss |], Theorem 3.1] showed 
that the expected number of leaves in a "pruned" 3x + 1 tree of depth d is 



We now briefly describe the search method to find the maximal value of a for a 
given level I, the largest number of ones allowed in a critical path in a certificate. 
We initialize the search with a taken to be the maximal value at level 1 — 1, and 
with the given certificate. Given a current test value of a, a rational number, we 
proceed with the algorithm above, looking for a certificate at level I. If one exists, 
we determine the critical path in the certificate having the largest ones-ratio, which 
gave a value a', which is our current champion. Then we take as a new test value 
a = a' + YQjjQQ and search it to level I for a certificate. If we get a new certificate 
by level I, we continue. Otherwise, if some tree remains unclosed at level I, we 
halt the computation, and the current value of a' is maximal at level This 
search procedure works because the maximal value a must be a rational number 
with denominator no larger than the maximal depth k < — , which in our search is 
always less than 100. All such values fall inside the Farey sequence of order 100, 
and the members of this Farey sequence differ by more than 10 q 00 ■ 

The straightforward method of computation to determine the critical path for 
a given 3x + 1 tree is to compute the entire tree. If the maximal level I to be 
searched is known in advance (as it is in the incremental search algorithm above) 
then we can get a speedup by pruning all branches of the tree at the point where 
they cannot contribute any critical path, i.e. even if extended to the full allowed 
depth k — L^J with all edge labels 1 they never have ones-ratio exceeding a. 
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To conclude the computational results, we give statistics on the structure of 
the certificates found, as described by the number of open vectors remaining at a 
given level I' < I of the certificate search; this value counts the number of prefixes 
in the certificate having length at least I'. We present data on these values for 
different values of a. This data is plotted on a logarithmic scale in Figure II When 




2 4 6 8 10 12 14 16 18 20 22 24 26 28 



level Z 

Figure 3. Number of open vectors in C by depth, for various 
values of a 

extrapolated, this data makes some indication of what depth of tree one is likely 
to have to search to get a certificate for larger values of a than we obtained. The 
rate of increase of a must slow down as I increases, so we view this extrapolation as 
indicating a lower bound on the depth to which one will have to search. We note, 
however, that the data in Figure || exhibits some irregular jumps the value of a, 
including a strange non-increase in the value of a from level Z = 5 to Z = 6. 

5. Discussion 

(1) How large a lower bound can the method of this paper produce? 
The stochastic models for 3a; + 1 trees (branching random walk models) studied 
in Lagarias and Weiss H predict that with probability one a random tree grown 
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to a large depth k — > oo will contain a path V having ^(V) > "/bp — e, for any 
fixed e > 0. However the method of proof given above requires the lower bound 
7o to be achieved for all trees at some given depth [], and not just for "almost 
all" trees of that depth. In terms of the stochastic models of || and j|, we are 
sampling an exponential number of trees (namely independently drawing 3 fe trees 
at depth k) out of the double-exponential number of possible trees in the branching 
random walk model at this depth, and demanding that they all have 7(7') > 7o- 
In our deterministic situation, if this event ever occurs at some level k , it yields 
a rigorous proof of the lower bound. In the analogous stochastic model, we should 
expect this cutoff event to occur at some finite level k with probability one whenever 
the set of "exceptional trees" with 7 < 70 has an exponentially small probability, 
less than 3~ fe , for all sufficiently large k. In || a large deviations analysis based 
on exactly this idea produced the constant 7sp. We conclude that the stochastic 
model supports the heuristic that this proof method should be expected to work for 
an Y 7o < Ibp, provided that we can search trees to a sufficient depth, depending 
on 70. 

To gain more confidence in this heuristic, it would be interesting to make the 
stochastic model analysis of || more precise, and to obtain a prediction for each 70, 
of what depth k one should expect to have to go to before an independently gen- 
erated sample of 3 fe trees of that depth k for the branching random walk, contains 
no maximal path with j(V) < 70 with probability at least ^. 

(2) Can a lower bound on the ones-ratio of | be attained? 

It remains surprising (to us) that a search to depth 60 of all trees was insufficient 
to produce a rigorous proof of the bound 

70 > r^TJ ~ 6.95212, 
log 4/3 

a bound which is believed to apply to the vast majority of initial values n. The 
statistical evidence provided by the cutoff values for increasing values of the ones 
ratio from 0.41 through 0.48, when extrapolated, seem to indicate that a search to 

3 We actually require only that each such tree contain some subtree containing the root vertex, 
of possibly smaller depth, having a path which achieves the given bound 7 > 70 
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max-weight I of at least 38, which corresponds to trees of depth 76, will be needed. 
We experimented on growing a few trees starting from the certificate with ones- 
ratio 0.48 to see how deep they had to become before a path with ones ratio of 0.50 
or greater appeared; these seemed to require depth around 68. If this represents 
the "knee" of the curve for 0.50 in Figure ||, then perhaps the depth necessary for 
a certificate might be smaller, say 72. In any case this would require a very large 
computation. 

With regard to extending the computations further, we note that the certificate 
search algorithm of §4 is well suited to parallel and distributed computation, since 
each open vector can be analyzed separately. Thus it should be possible to search 
to considerably greater depth using a network of machines. We leave a rigorous 
proof of the lower bound i as a challenge to future researchers. 
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